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We present a strategy for increasing thermoelectric performance by invoking ferroelectric distortion in a transition- 
metal oxide. By using an ab initio calculation approach with the Boltzmann transport equation, a large Seebeck co¬ 
efficient is calculated in rc-type BaTi0 3 . The polar structural distortion causes the peculiar dispersion of the lowest 
conduction band, which enhances the Seebeck coefficient along the polar direction. A microscopic mechanism of car¬ 
rier concentration dependence of the Seebeck coefficient is discussed in terms of the band velocity distribution in the 
Brillouin zone. 


1. Introduction 

Due to the recently growing energy issues, thermoelec¬ 
tric power has been attracting much attention since it can re¬ 
cover waste heat and be used as Peltier cooler without pro¬ 
ducing greenhouse gas emissions. As motivated by the ad¬ 
vances in the thermoelectric applications, it has been a key 
issue to improve the thermoelectric efficiency in material- 
science research field. The maximum thermoelectric perfor¬ 
mance at a temperature T is determined by a dimensionless 
figure of merit, ZT — S 2 ctT/k, which depends on the See¬ 
beck coefficient ( S ), the electrical conductivity (cr), and the 
thermal conductivity ( k ) including the electronic and lattice 
contributions. By assuming a parabolic electronic band and 
energy-independent scattering approximation, 1 S is described 
as S = (87rkg/ 3eh 2 )m*T(n/3n) 2 ^ 3 , where n is the carrier con¬ 
centration and m* is the effective mass of the carrier. The ther¬ 
mal conductivity is further decomposed into the electronic 
and phonon contributions: k-k s \ + K p h- Since the S (oc rT 1 ^) 
and cr (oc n) show the trade-off dependency on n, it is not easy 
to optimize the thermoelectric performance in the range of n 
limited for the practical use. Besides, m* may cause another 
conflict, such that larger effective mass produces higher S but 
lower cr. Therefore, designing high thermoelectric materials 
is to find the optimized balance of S, cr, and k . 23 The conven¬ 
tional strategy for enhancing ZT is to employ heavy elements 
and narrow-gap semiconductors as reducing k. For example, 
Bi 2 Te 3 —now popular as a representative topological insu¬ 
lator 4 — shows ZT > 1 at appropriate temperatures. 5 These 
materials are in practical use, however toxic to humans. 

After a striking finding of large Seebeck coefficient in 
Na A CoC> 2 , 6 transition-metal oxides have emerged as an¬ 
other category of thermoelectrics. Electron-doped SrTi0 3 and 
KTa0 3 have been found to have the large Seebeck coeffi¬ 
cients S , 7 ' 8 which are comparable to that of BhTei, as also 
having the high (metallic) conductivity cr. Several density- 
functional-theory (DFT) calculations have been performed 
for the thermoelectric transition-metal oxides, as they have 
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mostly concluded that the narrow t 2 g bands with the large m* 
are responsible for the large Seebeck coefficient. 9 - 10 Later on, 
Kuroki el al. have pointed out that a particular pudding-mold 
dig band may enhance both S and cr in dumbbell-type struc¬ 
ture in Na v Co02 and Q 1 AIO 2 . 11 ' 12 As for SrTi0 3 and KTa0 3 , 
they claimed that the multiplicity of three-fold t 2 g state plays 
a role in enhancement of S and cr. When the electron is doped 
into the empty t 2 g bands, the multiple bands can keep low 
Fermi level, which in turn results in a large Seebeck coeffi¬ 
cient, while the doped electrons increase the conductivity. 13 
Recently, Shirai et al. casted another interpretation on the 
large thermopower in SrTi0 3 as following. 14 It is often said 
that the t 2 g components of the d electrons form dispersion¬ 
less and three-fold degenerate bands. However, the lift of the 
degeneracy as leaving from the T point should not be over¬ 
looked. As far as the electron doping is concerned, only the 
lowest band is relevant. More importantly, even a single band 
can have electrons of quite different masses in different re¬ 
gions of the Brillouin zone. The fight electron leads to large cr, 
whereas the heavy mass results in large S, so as to increase the 
power factor S 2 cr. If this idea works, it is more critical to have 
strongly anisotropic m* in the single band rather than multi¬ 
ply degenerate bands for enhancement of the thermopower 
in oxides. In this context, we focus on the ts g bandstructure 
in BaTi0 3 , in which the ferroelectric crystal distortion is ex¬ 
pected to enhance anisotropy of the m* and to enhance the 
thermoelectric properties. In fact, an early experimental find¬ 
ing of the anisotropic Seebeck coefficient in BaTi0 3 was re¬ 
ported in 1967, 15 whereas the accurate theoretical approaches 
to understand the electronic property had not been available 
yet in that time. In the present study, we aim to confirm the 
effect of the polar distortion on the f 3g band and the resulting 
thermoelectric properties in BaTi0 3 and then we discuss the 
microscopic mechanism of enhancement of S and cr by means 
of first-principles DFT calculations. 

2. Bandstructure and Wannier Interpolation 

Perovskite BaTi0 3 is a popular ferroelectric oxide, which 
undergoes a ferroelectric transition at 393K. To obtain the 
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Fig. 1. Tetragonal unit cell of BaTi 03 in the (a) paraelectric (PE) structure 
and (b) ferroelectric (FE) structure. Apical and side oxygens are denoted as 
O ap and 0 s , respectively. Directions of the polar atomic distortions and the 
polarization are shown by block arrows. The polar displacements are exag¬ 
gerated by the factor of two for clarity. 



Fig. 2. Polar ionic displacement, defined as a change in bond length be¬ 
tween Ti and O ap ions in BaTiC> 3 , as a function of the electron doping rate 
x. 


electronic structure of ferroelectric BaTi 03 , we performed 
DFT calculations by using the VASP code 16 with the general¬ 
ized gradient approximation Perdew-Burke-Ernzerhof (GGA- 
PBE) potential. 17 The tetragonal lattice constants a=3.991A 
and c-4.035A were taken from experimental result at room 
temperature. 18 The internal atomic coordinates were opti¬ 
mized until forces acting on atoms were less than 1x10 3 
eV/A both for the paraelectric (centrosymmetric) structure 
and ferroelectric (polar) structure as shown in Fig.l. A (12, 
12, 12) Appoint mesh was used for the Brillouin zone inte¬ 
gration. Although the real crystal structure at the paraelectric 
phase (T > 393K) is known as cubic, we used the same tetrag¬ 
onal lattice for both paraelectric and ferroelectric structures to 
discuss the effect purely caused by the polar ionic displace¬ 
ment. The ferroelectric polarization in non-doped BaTi 03 
was calculated by Berry phase approach 19 as P DFT -26/uC/cm 2 
along the c axis, which is in good agreement with the exper¬ 
imental value of F >exp =27 / uC/cm 220 and the previously calcu¬ 
lated values of / jDFI -22-29pC/crrr. 21 ■ 22 The polarization can 
be decomposed into the ionic and electronic contributions: 


p DFT = P" m + p elec , where the ionic contribution was calcu¬ 
lated by a point-charge model as assuming the nominal ionic 
valence as +2, +4, -2 for Ba, Ti, O ions, respectively. The 
calculated P lon =14 /./C/cm 2 , which is almost half of the total 
polarization, originates from the Ti and O ionic displacement 
due to the strong hybridization between occupied O -p and un¬ 
occupied Ti-c/° state as called d°-ness mechanism. 23 Figure 2 
shows the polar ionic displacement as a function of electron 
doping. When electron is doped in BaTi 03 , the Fermi energy 
moves in the unoccupied Ti-t 2 g bands. It screens the electric 
potential of ionic charge and spoils the ferroelectricity. Nev¬ 
ertheless, the polar ionic displacement (i.e. the ionic contri¬ 
bution to P) is rather robust so as to stay almost 90% even 
when 0.05 e of electron is doped in the Ti-/f° state, as con¬ 
sistent with a previous theoretical study in Ref.24. This result 
ensures that the crystal structure keeps the polar distortion un¬ 
der some amount of electron doping (e.g. x < 0.1), which is 
sufficient for the following discussions on the electron-doped 
BaTi 03 in the ferroelectric structure. Technically speaking, 
the cubic-tetragonal (centrosymmetric-polar) transition tem¬ 
perature may decrease by electron doping although the polar 
ionic distortion will persist below the transition temperature. 
It was reported that the transition temperature is kept over 
room temperature until 0.022 e/u.c. doping, while higher dop¬ 
ing may lead to lower transition temperature. 25 For the ther¬ 
moelectric application of BaTi 03 , tuning of the polar property 
under electron doping will be required. For example, both the 
polar distortion and the transition temperature can be largely 
increased by epitaxial strain in BaTi 03 thin film. 26 - 27 

Figure 3 shows the bandstructure of BaTi 03 in the para¬ 
electric and ferroelectric structure. They consists of nine O-p 
bands below the Fermi energy, three Ti-fz# and two e g bands 
above the Fermi energy. Three T'i-? 2 ^ bands lying at the con¬ 
duction band bottom are well isolated from the other bands. 
Due to the tetragonal symmetry of the crystal structure, the t 2 g 
bands are split into (yz, zx) and xy bands at the T point. The 
former can be further split along the T-X line, where the four¬ 
fold rotation symmetry is broken. The symmetry-lowering ef¬ 
fect is significantly enhanced in the ferroelectric structure. It 
is intriguing that the (yz, zx ) bands are the conduction bottom 
in the paraelectric structure, whereas the xy band comes to the 
bottom in the ferroelectric structure. This is caused by strong 
Ti-O ap hybridization, which pushes up the anti-bonding (yz, 
zx) bands in the ferroelectric structure. In Fig.3 (d), it is ob¬ 
served that the xy and zx bands (which belong to the same 
irreducible representation A 2 ) repel each other along the T-X 
line as resulting in the particularly non-harmonic band shape. 
The lowest t 2 g band is responsible for the thermoelectric prop¬ 
erties since the electron will be doped into this band. For the 
further analysis of the thermoelectric property, the t 2 g bands 
were interpolated by means of maximally localized Wannier 
functions (MLWF) by using the Wannier9® code. 28,29 

Figure 4 shows one example of calculated MLWF for T\-d zx 
state. As related. Table I shows the hopping integral between 
Ti-d zx and O-p states as calculated via the off-diagonal Hamil- 
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Fig. 3. Bandstructure of BaTi 03 in (a) paraelectric and (b) ferroelectric 
phase, (c) and (d) are their blow-up versions. The Wannier-interpolated Ti- 
bands are superposed on the DFT bandstructure, as highlighted by (red) 
solid lines in (a) and (b). 



Fig. 4. An isosurface of MLWF for Ti- d zx orbital state hybridizing with 
the surrounding O-p states in the ferroelectric structure. The hopping integral 
between the d zx and the side (s) and apical (ap) oxygen p states are shown by 
arrows. 


directly related to the t/°-ness microscopic origin of ferroelec¬ 
tric distortion, by which the energy gap is increased to lower 
the band energy and to stabilize the ferroelectric structure. 23 
The anisotropic hopping integral slightly modulates the band- 
structure but has a great impact on the thermoelectric prop¬ 
erty. 


3. Thermoelectric Property 


The standard Boltzmann’s transport approach was adopted 
in the present study to calculate the Seebeck coefficients. 30 In 
this approach, the Seebeck coefficient is given by 


S(p,T) 


1 KiQi,T) 
eT K 0 (p, T) ’ 


(1) 


where e is electron charge, T is temperature, Kq and K\ are 
given by 



f 


b 

F 


/IT 

PE 

1.03 

0.98 

0.98 

2.00 

2.02 

2.02 

FE 

1.03 

1.34 

0.72 

2.00 

1.89 

2.15 


Table I. Calculated hopping integrals t (eV) and the bond length / (A) be¬ 
tween Ti-dtyZ'Zx) and side (s) and apical (ap) O-p orbital states in the paraelec¬ 
tric (PE) and ferroelectric (FE) structures. 


tonian element for the MLWF basis together with the Ti-O 
bond length. In the ferroelectric structure, a Ti ion is slightly 
displaced toward one of the apical oxygens by 0.13 A, which 
considerably enhances the hopping integral. Therefore, the 
MLWF of Ti-d zx has larger weight at the top apical oxygen 
site than the bottom apical oxygen site. As shown in Table I, 
the calculated hopping with closer O-p state is almost double 
of the other (f“ p >> L‘ p ). The anisotropic enhancement of t is 


Up, T) 



TVnkVnk 


df(s,p, T ) 
ds 


0 Snk-PY ( 2 ) 


for ; =0,1. Here, s nk is the band energy, v nk is the band ve¬ 
locity, r is the quasiparticle lifetime assumed as constant, /(e) 
is the Fermi distribution function, and p is the chemical po¬ 
tential. In order to calculate the accurate band velocity (with 
band index n) 


Vnk = T 


1 ds nk 


H dk 


(3) 


Wannier function approach implemented in BoltzWann 31 
code was employed. The Appoint mesh was increased up to 
(160, 160, 160) to reach the convergence. The electrical con¬ 
ductivity cr is also calculated as 


ir(p, T) 


eh r 

— ,L 


de(- 


df(s,n , T ) 
ds 


>E 

n, k 


v lk 6 ^- 


(4) 
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Fig. 5. (a) Calculated Seebeck coefficient for x=0.03 plotted as a function 

of temperature for ferroelectric (FE) and paraelectric (PE) crystal structures, 
(b) Calculated Seebeck coefficient and (c) power factor (PF) (devided by t) 
as a function of the doping rate x at T =300 K. The diagonal matrix elements 
zz (the polar axis) and jo: are plotted by solid and dotted lines, where broad 
and narrow lines correspond to FE and PE structures, respectively. 


In Fig.5 (a), we show the calculated Seebeck coefficients as 
a function of temperature T at a fixed doping rate x=0.03. In 
this calculation, the rigid band approximation was assumed. 
In the paraelectric structure, S X x and S zz show the similar val¬ 
ues. S xx reaches -122 pV/K at T=300K and x=0.03, being 
consistent with the previously calculated value for SrTi 03 as 
Sxx = -87pV/K at T =300K and x=0.05. 13 However, as it was 
discussed in Ref. 13, the calculated S value for SrTiCL is al¬ 
most half-reduced from the experimental result. This is due 
to the strong correlation effect of the 3d state, which band 
width is usually overestimated in DFT calculations. As for 
the further correspondence to the real value of S , we should 


note that the S value predicted by DFT calculations is under¬ 
estimated roughly by the factor of two. For the ferroelectric 
case, the calculated Seebeck coefficient shows the remarkable 
anisotropy, S zz » S xx . At T=300K and x=0.03, S zz is cal¬ 
culated as -230pV/K, which is almost four times larger than 
S xx — -74 pV/K. This value may be comparable to the ex¬ 
perimentally measured S =s300pV/K in epitaxial La-doped 
BaTiCL. 32 Figure 5 (b) shows the calculated Seebeck coeffi¬ 
cient as a function of the doping rate x. The .S' „ coefficient be¬ 
haves anomalously as increasing x, as making a large “bump” 
at x=0.03. This is caused by the anisotropic property of the 
Fermi surfaces and the effective mass due to the ferroelec¬ 
tric ionic displacement, as explained in later section. It results 
in a remarkable crossover of power factor (S 2 <x) between the 
paraelectric and ferroelectric phases at x ~0.05 as shown in 
Fig. 5 (c). 

4. Effect of Ferroelectric Distortion 

Here, we discuss the microscopic origin of the enhance¬ 
ment and the anisotropic behavior of Seebeck coefficient S 
in ferroelectric BaTiC> 3 . We note that the enhancement of S zz 
originates from a single-band contribution. In the energy re¬ 
gion of 0 < p < 0.5 eV, a main portion of the electron density 
of states originates from the xy band at the conduction bot¬ 
tom, whereas a contribution from the (yz, zx) bands is less 
than 5% as being irrelevant to the trend of S(j-i). This makes 
a clear contrast to the cubic SrTiC >3 case, where the three t 2s , 
bands have the equal contribution to S(p). 11 

As already discussed by Kuroki et al. , the magnitude of 
Seebeck coefficient is mainly determined by the energy de¬ 
pendence of the electron velocity v„k- 13 By making a rough 
approximation, Kq and K\ in Eq.(2) are described as 

Kq ~ ^ ^ (^high + Vlow“), ~ k H T ^ ^ (Vhigh~ — H ow ), (5) 

where the summation is over the states in the range of \sk - 
p| < kuT and Vhi g h and vi ow are typical velocities for the states 
above and below p, respectively. Therefore, Vhigh » vi ow situ¬ 
ation leads to a large Seebeck coefficients. 

Figure 6 (a) and (b) show the calculated v 2 ^ superimposed 
on the sections of iso-energy surfaces e n k with several p val¬ 
ues. As increasing /r, the ellipsoidal energy isosurface is elon¬ 
gated along the k z direction for /./ < 0.2 eV, and then the lobes 
grow up along the k x and k y direction for /j > 0.2 eV. By tak¬ 
ing a close look at Fig. 6 (a) and (b), the ellipsoidal energy 
surface has some amount of v 2 but almost zero v 2 . When /j > 
0.2 eV, the energy surface starts to have finite v 2 value. There¬ 
fore, the expectation value of (v?)(//) is shifted rightward by 
0.2 eV with respect to (v 2 )(yu) as shown in Fig. 6 (c). For /r ~ 
0.2 eV, the situation V; lgh » vl ow is realized so as to signifi¬ 
cantly increase S zz . This can explain the “bump” behavior of 
5^(x) in Fig. 5 (b) at x=0.03 (the relation between /r and x is 
shown in Fig. 6 (d)). 

The “bump” behavior of S (x) has been recently reported 
with the experimental observation in SnTe. 33 SnTe is now 
known as a topological crystalline insulator, where the Sn-p 
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Fig. 6. Contour plot of s(k)=ii surfaces with several /i values (origin of the 
chemical potential /i is taken as the bottom of the conduction bands) and the 
squared band velocity plotted by color for (a) v\ and (b) v? of the lowest con¬ 
duction band (xy band) in ferroelectric BaTi 03 . The three-dimensional en¬ 
ergy surfaces are shown in the right panel (i)-(iii). (c) Squared band velocity 
(v 2 x ) and <v 2 ) integrated in the k -space as <v,y) 2 (/r) = (1/VO Z„,k v in kvjnkSip - 
£ n k) for the lowest conduction band, (d) Chemical potential fi as a function 
of the doping rate x at T =300K. 


and 'Yc-p orbital characters are inverted between the valence- 
band top and conduction-band bottom at the R point. 34 35 
When a hole carrier is doped, this band inversion causes the 
warping in the Fermi surfaces of the valence band, which 
in turn leads to the enhancement of S (x). 10,36 In short, 
in SnTe and the related system, the strong spin-orbit cou¬ 
pling mixing the orbital character between narrow gap causes 
non-harmonic Fermi surfaces, which is responsible for the 
high thermoelectric property. This is somehow related to the 


present case of BaTiC> 3 , where the polar crystal distortion 
asymmetries the Fermi surfaces and results in the sizable en¬ 
hancement of S (x) along the polar direction. By using this 
unique property, it may be possible to control the thermoelec¬ 
tric performance by the piezo-electric effect. When one im¬ 
poses uniaxial strain on BaTiC> 3 , both the polar axis and the 
thermoelectrically favored direction may be manipulated as 
one wishes. This effect needs to be confirmed by further the¬ 
oretical and experimental studies. 

5. Conclusion 

In summary, within the Boltzmann transport theory we 
calculated an anisotropic Seebeck coefficient in ferroelectric 
BaTiCh. The polar displacements of Ti and O ions cause the 
t 2 g band splitting into the (yz, z.x) and the xy states: the lat¬ 
ter band in turn results in the strong anisotropy and enhance¬ 
ment of the Seebeck coefficients. Such a structural distortion 
in perovskite oxides can provide a way to asymmetrize the 
constant energy surfaces and the unusual energy dependence 
of the band velocity. This finding may pave a way for explor¬ 
ing novel thermoelectrics in various ferroelectric oxides. For 
example, PbTiC >3 may cause a larger Seebeck coefficient as 
having the huge electric polarization. 
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